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Abstract 

Nuclear level densities for ^^'^'^Fe have been extracted from the primary 7-ray spectra using 
(^He,^He'7) and (^He,a7) reactions. Nuclear thermodynamic properties for ^^Fe and ^''Fe are 
investigated using the experimental level densities. These properties include entropy, Helmholtz 
free energy, caloric curves, chemical potential, and heat capacity. In particular, the breaking of 
Cooper pairs and single-quasiparticle entropy are discussed and shown to be important concepts for 
describing nuclear level density. Microscopic model calculations are performed for level densities 
of ^^'^''Fe. The experimental and calculated level densities are compared. The average number of 
broken Cooper pairs and the parity distribution are extracted as a function of excitation energy 
for ^^'^''Fe from the model calculations. 
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I. INTRODUCTION 



Nuclear thermodynamics has attracted considerable attention in recent years. Several 
temperature-dependent nuclear properties such as nuclear shapes, giant dipole resonance 
widths and their fluctuation properties have been investigated in the literature l|. In this 
context one of the most interesting topics is that of phase transitions in atomic nuclei. 
Different types of phase transitions have been discussed in nuclear physics. One is the flrst 
order phase transition in the multifragmentation of nuclei which occurs at high temperatures 
[2I. Negative heat capacities observed in multifragmentation of nuclei have been obtained 
Tom energy fluctuations and interpreted as an indicator of a flrst order phase transition 

0,Q. 

The second type of phase transition in atomic nuclei is the transition from a phase with 
strong pairing correlations to a phase with weak pairing correlations 4|. The onset of a 
discontinuity in thermodynamic variables and the evolution of zeros of the canonical [5|, 16| 
and grand-canonical partition functions in the complex plane have been discussed in 
terms of pairing transitions. 

Recently [tI, thermal properties for ^^Fe have been calculated using both static-path 
approximation (SPA) and SPA plus random-phase approximation (RPA). These calculations 
show that the increase of the moment of inertia with increasing temperature is correlated 
with the suppression of pairing correlations. 

Furthermore, structures in the heat-capacity curve related to the quenching of pairing cor- 
relations have been obtained within relativistic mean fleld theory Jsl, the flnite-temperature 
iartree-Fock-Bogoliubov theory j^, and the shell- model Monte Carlo (SMMC) approach 
lo| . An S-shaped structure in the heat capacity curve derived from the level densities of 
low-spin states has also been observed experimentally ll|] and interpreted as a signature of 
a pairing phase transition. 

As a part of an ongoing effort, in the present work we extract several thermodynamic 
properties of ^^Fe and ^^Fe isotopes starting from nuclear level densities. Details of the 
experiment, analysis tools, the Oslo method, and experimental level densities are presented 
in Sect. [Tll Thermodynamic concepts are discussed in Sect. IIIII and comparison with a 
combinatorial model is shown in Sect. IIV[ Finally our results are summarized in Sect. |Vl 
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II. EXPERIMENTAL DETAILS AND DATA ANALYSIS 



The self-supporting ^''Fe target was bombarded by a 2-nA beam of 45-MeV ^He particles 
from the Oslo Cyclotron Laboratory at the University of Oslo. The target was isotopically 
enriched to 94.7%, and had a thickness of 3.4 mg/cm^. The outgoing charged particles were 
recorded by eight Si AE-E telescopes, which were collimated and placed 5 cm away from the 
target at a ring located at 45° with respect to the beam direction. The particle telescopes 
covered 0.3% of the total solid angle. The thicknesses of the front and the end detectors 
were 140 and 3000 /^m, respectively, and the particle energy resolution was ~ 0.3 MeV over 
the entire spectrum. The reaction 7 rays were detected by 28 collimated 5" x 5" Nal(Tl) 
detectors with a total efficiency of ^ 15% of Atc, and with 6% energy resolution at 1.3 MeV. 
In order to monitor the selectivity and populated spin distribution of the reactions, one 60% 
Ge(HP) detector was used. 

The excitation energy of the final nucleus is determined from the known Q value and 
the reaction kinematics using particle-7 coincidence information. A total 7-ray spectrum is 
obtained for each excitation energy region. These 7-ray spectra are then unfolded using a 
Compton subtraction method 12|]. A two-dimensional primary 7-ray matrix is extracted by 



applying a subtraction method to the unfolded 7-ray spectra |13| . Basic assumptions and 



details of the subtraction method are given in Ref. 

The primary 7-ray matrix is factorized using the Brink- Axel hypothesis ,[l5|, according 
to which the probability of emitting a 7 ray from an excited state is proportional to the 
7-ray transmission coefficient T{E^) and the level density at the final energy p{E — E^). 
This factorization is determined by a least method without assuming any functional form 



for the level density and the 7-ray transmission coefficient 16|- However, this method does 
not provide a unique solution. This factorization is invariant under the transformation [l^ 



p{E-E^) = Aexp{a{E - E^))p{E - E^) 

f{E^) = Bexp{aE^)T{E^), (1) 

where A, B, and a are the free parameters of the transformation. Therefore it is very 
important to determine accurately the free parameters A, B, and a in order to find the 
physical solution. The parameters A and a are determined from the normalization of the 
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level density to the discrete levels at low excitation energies and to the density of the neutron 
resonances at the neutron binding energy Bn- The parameter B is then determined using 
the average total radiative width of neutron resonances jl?! . 

Since there are no neutron resonance data for the ^^Fe compound nucleus a different 
procedure for the normalization of the level density is performed. We use ^^Fe as a basis, 
since both average neutron resonance spacings D and total radiative widths of neutron 



resonances are well known [l8|]. The procedure for extractin g th e total level density p(-Bn) 
from the average resonance spacing D is described in Ref. 16|. Then we apply the von 



Egidy and Bucurescu parameterization of the back-shifted Fermi gas formula 
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exp(2v^aZ7) 

^""^^ = ^2v/2aV4t/5/v' (2) 

where a is the level density parameter, and the intrisic excitation energy U = E — Ei is 
determined by the backshift parameter Ei. In this global description of the level density. 



shell corrections are taken into account in the estimation of the a and Ei parameters 19|. 
By fitting Pbsfg to p{Bn) in ^^Fe, we determine the normalization parameter rj = 0.852 
of Eq. ([2]). This value is then adopted for ^^Fe together with the prescribed a and Ei 
parameters. The parameters used are listed in Table [H 

Figure [1] shows the level densities of ^^•^'^Fe from the ground state up to ~ Bn — 1 MeV. 
They are normalized to discrete levels at low excitation energies (jagged lines) and to the 
level densities at the neutron binding energy Bn (triangles). The arrows indicate the fitting 
regions used. As described above, the triangles are determined from Eq. ([2]) with t] = 0.852 
for ^^Fe and from the observed D value for ^^Fe. 

The level density obtained with the Oslo method follows well that obtained by counting 
of discrete levels up to around 5 and 2.5 MeV of excitation energies in ^^'^''Fe, respectively. 
Both nuclei show an abrupt increase in level density ai E = 2 — 3 MeV, indicating the first 
breaking of nucleon Cooper pairs. 

Figure [1] also includes level densities from particle evaporation data (open circles) mea- 
sured withjthe 55Mn(d,n)56Fe reaction and with the 5SFe(3He,a)^^Fe and ^^Co{d, af^Fe 



reactions [2l|]. The method determines the slope of the level densities, but not the absolute 
normalization constant. Thus, these data have been scaled to match pbsfg (solid smooth 
curves) . The slopes of the level densities from evaporation studies fit very well with the em- 
ployed Pbsfg- This gives support to the level density parameters a extracted as prescribed 
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FIG. 1: (a) Experimental level density of Fe (full circles). The jagged solid lines represent 
the level density obtained from counting of discrete levels 



23( 1 . The smooth solid curve is the 
Q. The data 



renormalized level density parametrization according to von Egidy and Bucurescu 
points between the arrows at low and high excitations are used for the normalization. The triangle 
at Bn is determined from pnsFG with ry = 0.852. The level density obtained from the ^^Mn(d,n)^^Fe 
reaction [2(i] (open circles) is scaled with a factor of 1.3. (b) Experimental level density of ^^Fe 
(full circles). The triangle at Bn is determined from neutron resonance spacings. The level density 



obtained from the Fe( He,a) Fe reaction 2J| (open circles) is scaled with a factor of 2.6 



in Ref. m and listed in Table [B 

From Eq. ([1]) we see that the a parameter determines the slope of the transmission 



coefficient, and thereby also the radiative strength function (RSF) f{E^) = T 



iE^)/{2^El 



assuming only dipole radiation. The RSFs have been published previously [20|, , but with 
other normalization procedures. Thus, for completeness, we show in Fig. [2] the renormalized 
RSFs. In the ^^Fe case, the B parameter of Eg . Q was determined from the average 
radiative width (F^) of neutron resonances at [l8|. For ^^Fe, where no (F^) exists, we 
scale the total RSF to match ^^Fe. The parameters are listed in Table [H 
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FIG. 2: Radiative strength functions for ^^Fe and ^'^Fe. 
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TABLE I: Back-shifted Fermi gas level density parameters and normalization con- 
stants. 



Isotope 


Bn a 


El D a{Bn) 


V 








(McV) (McV-^ 


) (McV) (k(A') 




(McV-i) 


(mcV) 


56pg 


11.197 6.196 


0.942 2.6(6)" 4.049 


0.852« 


3600(800)'^ 


2300(1200)" 


57Fg 


7.646 6.581 


-0.523 22.0(17) 3.834 


0.852 


1380(170) 


900(470) 



" Based on systematics, see text. 



III. THERMODYNAMIC QUANTITIES 



Depending on the system under study, one can chose among different kinds of statistical 
ensembles in order to derive thermodynamic quantities. The thermodynamic quantities 
derived within different ensembles give the same results in the thermodynamic limit. On 
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the other hand, the choice of a specific ensemble may change results significantly for small 
systems. For example, the caloric curves derived within the microcanonical and canonical 
ensembles coincide for large systems; but the two caloric curves depart from each other for 



small systems 



24, 



25|. 



The microcanonical ensemble is commonly accepted as the appropriate ensemble to use 
in investigating atomic nuclei, as the nuclear force has a short range and the nucleus does 
not share its excitation energy with its surroundings. However some thermodynamic quan- 
tities such as temperatures and heat capacities may have large fiuctuations and negative 
values when derived within the microcanonical ensemble. On the other hand, the canonical 
ensemble averages too much over structural changes of the system. Therefore, it is diffucult 
to chose an appropriate ensemble for a small system. Here, we use both microcanonical and 
canonical ensembles to study the thermodynamic properties of the system. 



A. Microcanonical Ensemble 



The microcanonical entropy is closely related to the level density of the system at a given 
excitation energy. Several thermodynamic properties of the atomic nucleus can be derived 
from the entropy. The entropy is defined as the natural logarithm of the multiplicity Q of 
accessible states within the energy interval E and E + 6E: 

S{E) = kBlnQ{E), (3) 

where kB is the Boltzmann's constant. Here we set to unity so that the entropy becomes 
dimensionless and thus the temperature T has the unit of MeV. One should note that the 
experimental level density is not the true multiplicity of states; i.e., it does not include the 
(2/ + 1) degeneracy of magnetic substates. Thus, in order to obtain the state density Q, one 
needs to know the spin distribution as a function of excitation energy. The spin distribution 
is usually assumed to be Gaussian with a mean of {21 + 1) = ^/2TTa, where a is the spin 
cut-off parameter which depends very weakly (a oc E^^^) on excitation energy. 

However, in this work we define a "pseudo" entropy based on the experimental level 
density, i.e., Q{E) = p{E)/pQ. The normahzation constant po is introduced and adjusted 
such that the third law of thermodynamics is fulfilled as S{T — 0) = 0. The ground 
states of even-even nuclei represent a well-ordered system with no thermal excitation and 
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FIG. 3: Upper panel: Deduced entropies for ^^'^^Fe. The solid line is the constant temperature 
least square fit to the data. Lower panel: Deduced entropy excess for the single particle calculated 
using Eq. dSJ. The entropy excess is A5 = 1.2(4). 

are characterized by zero entropy and temperature. Therefore the normahzation constant 
is set to po = 1 MeV~^ in order to obtain 5 = lni7 ~ in the ground state band region of 
^^Fe. The extracted po is also used for the ^^Fe nucleus. 

Figure [3] shows the entropies of ^^'^^Fe. The breaking of the first Cooper pair appears in 
the E = 2 — 3 MeV excitation energy region for both nuclei, indicating a slightly delayed 
breaking for ^^Fe due to the odd valence neutron. 

The entropies of ^^•^'^Fe appear fairly linear at high excitation energies, and the slope of 
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the entropy is related to the temperature of the system by 

T = {dS/dE)y\ (4) 

The entropies of ^^'^'''Fe are fit with a constant temperature model, shown with solid lines 
in Fig. [3l From this model, constant temperatures of T = 1.2(3) MeV and T = 1.0(3) MeV 
are found for ^^Fe and ^^Fe, respectively. These temperatures are interpreted as the critical 
temperatures Tc for the breaking of nucleon pairs. 

The entropy difference between the even-odd and even-even nuclei is interpreted as the 
entropy of a single quasiparticle (particle or hole), assuming that the entropy is an extensive 
(additive) quantity. The entropy carried by the valence neutron particle in our case can be 
estimated by 



AS = SC'^Fe) - SC'^Fe). (5) 

The lower panel of Fig. [3] shows the single particle entropy as a function of excitation energy. 
The fluctuations at the low excitation energies result from the lower pairing gap in the odd- 
mass system. At higher excitation energies above ~ 4 MeV, the entropy difference has a 
relatively constant value, revealing the statistical behavior of these isotopes, AS* = 1.2(4' . 
This value is less than the value obtained for the rare-earth isotopes which is AS* = 1.7 26| . 
This is expected because ^^Fe and ^''Fe isotopes are in the vicinity of closed shells, thus the 
entropy is less than that of the rare-earth isotopes. 

Assuming a constant AS* and a constant energy shift A between the two entropies, these 
two quantities are connected to the critical temperature Tc (26| by 

T. = ^A. (6) 

From this relation, we can calculate the entropy difference AS = A/Tc = 1.5/1.1 ~ 1.4, 
which is consistent with the present observations. Here, the gap parameter A is assumed 
to include both the pairing gap and the contribution from the distance between the Fermi 
surface and the closest neutron orbital. 



B. Canonical Ensemble 



In a canonical ensemble, a physical system is in thermal equilibrium with a heat reservoir 
at constant temperature and can exchange energy (but not particles) with the reservoir. 
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Analogous to the entropy in the microcanonical ensemble, the partition function is the 
starting point to obtain the thermodynamics of a system in the canonical ensemble. The 
partition function for a given temperature is determined by 



Z{T) = Y,PmEe-'^l^, (7) 



£=0 



where p{E) is the measured level density, and 5E is the energy bin used. The summation 
in Z{T) goes to infinity and our level densities extend to ~ Bn — 1 MeV. Therefore we 
extrapolate the experimental level densities using the BSFG model parameterized by von 



Egidy and Bucurescu 19|. The Helmholtz free energy can be calculated from the partition 
function by 

F(r) = -rinZ(T). (8) 

This equation gives the connection between statistical mechanics and thermodynamics in the 
canonical ensemble, as does the entropy in the microcanonical ensemble. Then the entropy 
S", the average excitation energy {E), the heat capacity Cy, and the chemical potential /z at 
a given temperature can be derived from F{T) by 

{E{T)) = F + TS, (10) 

MT) ^ f , (12) 

where n is the number of thermal particles. These thermal particles outside a core of Cooper 
pairs are responsible for the thermal properties of the nucleus at low excitations. At higher 
temperatures, pairing correlations are quenched, and the nucleus has a transition from a 
paired to unpaired phase. 

Figure m shows the Helmholtz free energy F, average excitation energy {E{T)), entropy 
5", and heat capacity Cy for ^^Fe (dashed lines) and ^^Fe (solid lines). The Helmholtz 
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FIG. 4: The Helmholtz free energy F, average excitation energy {E{T)), entropy S, and heat 
capacity Cy for ^^Fe (dashed lines) and ^^Fe (sohd lines), deduced in the canonical ensemble. 

free energy F and the average excitation energy {E(T)) behave smoothly as a function of 
temperature. At around T ~ 1.3 MeV, both isotopes are excited to energies comparable 
to their respective neutron binding energies. Average excitation energies for ^^Fe and ^^Fe 
coincide only at one point T ~ 2 MeV. Above this temperature the odd isotope has larger 
values as the temperature increases. The entropy S and heat capacity Cy are the first 
and second derivatives of F, respectively, and thus both reflect thermal variations. For 
temperatures below T = 1 MeV, the entropy difference for ^^Fe and ^^Fe reaches ~ 2, the 
entropy of ^''Fe being larger. The entropies have similar values for T = 1 — 1.5 MeV. Above 
T = 1.5 MeV, the two entropies start to diverge. Both heat capacities have an S shape 
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FIG. 5: The experimental Helmholtz free energies deduced in the canonical ensemble for ^"Fe and 
^^Fe. The critical temperature for the quenching of pairing correlations is indicated by T^, where 
the chemical potential is ~ 0. 

which is interpreted as a fingerprint for pairing transitions in nuclei. The structure in the 
heat capacity Cy for ^^Fe is more pronounced than the Cy for ^^Fe, which is the opposite 
of what SMMC calculations predict. The contribution to the heat capacity from collective 
excitations is negligible, and has no influence on the S shape [271] . 

Figure [5] shows the chemical potential which is defined as the amount of energy required 
to excite a nucleon from the underlying core of paired nucleons, while the entropy and volume 
held fixed. The typical energy cost for creating a quasiparticle is —A which is also equal to 
the chemical potential. The chemical potential can be written as 
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^=AN = 1 = ~^- ^'^^ 

thus giving F^yen = Fodd + A- The energy curve of ^''Fe can be interpreted as the free energy 
for an even-even system with one extra nucleon. 



C. Probability Density Functions 

An alternative way of studying the energy-temperature relation of a system is the prob- 
ability density function, which is the probability of the system having energy E for a given 
temperature, and is given by 

. 5(£)|Mz£/I). (14, 

In the thermodynamic limit P{E) is very sharp - similar to a 5 function. However, for small 
systems the multiplicity of states is much smaller, which makes the probability distribution 
broader. Figure [6] shows the probability density functions for ^^•^'^Fe nuclei for various 
temperatures. Below T ~ 1 MeV, the distribution is mainly based on the experimental 
level density. For increasing temperatures, the energy of the nucleus increases; therefore, 
the probability density function depends more and more on the extrapolated level density. 

Below T = 1.2 MeV, the distribution is slightly broader for ^^Fe due to collective exci- 
tations below the pairing gap in this nucleus. The effect of critical temperatures obtained 
with the constant temperature least-square fit can also be seen on the probability density 
function. At = 1.2 MeV and = 1.0 MeV for ^^'^^Fe, respectively, the probability 
distribution becomes broader where the breaking of nucleon pair process is strongest. 

IV. COMPARISON OF LEVEL DENSITIES WITH MICROSCOPIC CALCULA- 
TIONS 

In order to extract more information on the underlying nuclear structure resulting in 
the observed level density, we have performed microscopic calculations with the code MI- 
CRO [2^. The code is based on a model employing Bardeen-Cooper-Schrieffer (BCS) quasi- 
particles 2^, which are distributed on all possible proton and neutron configurations ac- 
cording to the given excitation energy E of the nucleus. The advantage of this code is a fast 
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FIG. 6: Observed probability density functions for ^^-^^pg fQy. various temperatures. 

algorithm that may include a large model space of single-particle states. Since level density 
is a gross property, the detailed knowledge of the many-particle matrix elements through 
large diagonalizing algorithms is not necessary. 

The single-particle energies Cgp are calculated with the Nilsson Hamiltonian for an axially 
deformed core with a quadrupole deformation parameter €2- The spin-orbit and centrifugal 
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TABLE II: Input parameters for the model calculations. 

Isotope £2 Att A^ ^lot ^^^0 ?ilJJvib Att A,y 

(MeV) (MeV) (MeV) (MeV) (MeV) (MeV) (MeV) 

56Fe 0.24 1.568 1.363 0.120 10.65 2.656 45.89 48.23 

"Fe 0.25 1.268 1.465^^ 0.120 10.72 2.656 45.66 48.40 
Taken from ^^Fe. 

parameters k and /i, together with the oscillator quantum energy 

between the harmonic oscillator shells, are also input to the code. Within the BCS model, 
the single-quasiparticle energies are defined by 



qp- ^,^sp-A)2 + A2, (15) 

where the Fermi level A is adjusted to reproduce the number of particles in the system and 
A is the pair-gap parameter. 

In the calculations we have adopted the Nilsson parameters k = 0.066 and fi = 0.32 
taken from Ref. ^Q\. The quadrupole deformation €2 was set to 0.24 for ^^Fe, and 0.25 for 
^^Fe jsij. The Nilsson levels used in the calculations for ^^Fe are shown in Fig. [TJ with the 
Fermi levels for the protons and neutrons. 

The rotational and vibrational terms, which are schematically added, contribute signifi- 
cantly to the total level density only in the lower excitation region. The rotational parameter 



Ay-ot was set to 0.12 MeV in order to reproduce the ground-state rotational band of ^^Fe 



23| 



The adopted pairing gap parameters A,r and A,^ were evaluated from even-odd mass differ- 



ences I34] according to Ref. [33| (see Table [ITl for a list of all parameters used) . 

The experimental and calculated level densities for ^^'^"^Fe are shown in Fig. [HI In general, 
there is a very good overall agreement between the measured level densities and the calcula- 
tions, even without reducing the pairing gap. This is in contradiction with what is usually 
thought using the canonical ensemble. However, excitation energy around 8 MeV may be 
too low to see the overall quenching of the pairing correlations demonstrated in Figs, (jl]) 
and (jS]). Microcanonical and canonical ensembles yield different results for small systems. It 
is, therefore, no surprise that one measures different thermodynamic quantities within these 
ensembles. 

In Fig. [HI especially for ^^Fe, it is gratifying how both the general functional form and 
the absolute magnitude of the level densities are reproduced. In the case of ^^Fe, there is 
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FIG. 7: Nilsson level scheme for ^^Fe. 

an overshoot in the calculated level density compared to the experimental data from about 
3.5 MeV excitation energy. This excess could be due to the neutron /5/2, pi/2 and (79/2 
orbitals coming into play too soon, that is, they are too close to the neutron Fermi level. 
Another possibility is that the adopted (constant) deformation is too large at these excitation 
energies, where it may be that the nuclear structure favors a more spherical shape. 

Figure [9] shows the average number of broken Cooper pairs (A^qp) as a function of excita- 
tion energy. Both neutron and proton pairs are taken into account, adding up to the total 
number of broken pairs. From Fig. [9], the pair-breaking process is seen to start a.t E 2.5 
MeV for both nuclei, in accordance with the values used for the pair-gap parameters (see 
Table HH). 
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FIG. 8: Calculated level densities (solid line) compared with the experimental ones (data points 
with error bars) for ^^'^'^Fe. 

A prominent step structure at (A'^qp) = 1 is seen clearly in ^^Fe for excitation energies 
between 2.5 < E < 5.5 MeV, and for ^''Fe in the energy region 4.0 < E < 7.5 MeV. This 
means that on average there is one Cooper pair broken, which could be either a proton 
or a neutron pair. For ^^Fe, the total (A^^qp) exhibits an increase for E > 5.5 MeV up to 
about 9 MeV, where there seems to be a saturation for (Aqp) ^ 2. This is also reflected 
in the proton and neutron contributions to (A"qp). A similar behavior can be seen in ^^Fe 
for excitation energies higher than ^ 8 MeV. We observe that the onset of this increase 
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FIG. 9: The average number of broken Cooper pairs as a function of excitation energy for ^"'^'Fe. 

is at higher excitation energies than for ^^Fe. A possible explanation of this is that the 
^^Fe valence neutron blocks the quasi-particles created during the pair-breaking process, 
suppressing the average number of quasi-particles even at high excitation energies. This is 
supported by the fact that the neutron contribution to the total {N^p) is lower or about the 
same as the proton contribution in ^^Fe, while in ^^Fe, the neutron contribution is in general 
higher than the proton contribution. 

The present model gives the opportunity of investigating the parity distribution as a 
function of excitation energy. For this purpose, we utilize the parity asymmetry defined in 
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Ref. y by 



a = , (16) 

P+ + P- 

which becomes —1 if only negative parity states are present, +1 if there are only positive 
parity states, and when both parities are equally represented. In literature, one also find 
the expression p-/p+ which relates to a by 



Pz. 

P+ 



1 — a 
1 + a' 



(17) 



The parity asymmetry a for ^^'^'^Fe is shown in Fig. [TOl 

From Fig. [101 we see that there is an excess of positive parity states in ^^Fe and negative 
parity states in ^''Fe. As the excitation energy increases, the parity asymmetry decreases for 
both nuclei, and shows an almost decay-like behavior, especially in the case of ^^Fe. Here, 
the parity asymmetry is close to zero at = 9 — 10 MeV. For ^^Fe, there is still a significant 
overshoot of positive-parity states in this energy region, giving an average parity asymmetry 
of a ^ 0.35 for 9.5 ~ E ~ 10 MeV. 

We have compared our results with the shell-model Monte Carlo results of Y. Alhassid 



et al. 
et al. 



35|, and with recent macroscopic-microscopic calculations performed by D. Mocelj 
361]. In Fig. 2 of Ref. 36|], the ratio p^/p+ is shown for ^^Fe, indicating a value of 



P-/p+ — 0.1 at 10 MeV excitation energy. From Fig. 4 in Ref. [35|], the ratio p^/p+ — 0.2 
for E = 10 MeV. Using an average parity asymmetry of a ~ 0.35 as determined in the 
previous section, we obtain p-/p+ ~ 0.5, which implies a considerable number of additional 
positive-parity states than predicted by the macroscopic-microscopic model and the shell- 
model Monte Carlo approach. For ^''Fe with a ~ 0.1 for 9.5 ~ E' ~ 10 MeV, we find 
P-/p^ ~ 0.8. However, one should note that in our calculations, the parity asymmetry 
strongly depends on the position of the neutron gg/2 orbital relative to the Fermi level (see 

Fig. ED- 



V. SUMMARY 

Nuclear level densities for ^^'^''Fe are renormalized using the new level density parameter- 
ization suggested by von Egidy and Bucurescu. The level densities obtained with the Oslo 
method agree well with those obtained from other experiments. The experimental level den- 
sities are used to extract thermodynamic quantities. The entropies for ^^•^'^Fe obtained in the 
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FIG. 10: The parity asymmetry as function of excitation energy for °°'°'Fe. 

microcanonical ensemble reveal step structures indicating the breaking of nucleon Cooper 
pairs. The entropy carried by the single neutron at higher excitation energies (4 MeV < E < 
7 MeV) is estimated to be AS — 1.2(4)A;b which is smaller than that of the rare-earth iso- 
topes. 

Assuming a constant AS, a critical temperature for the depairing process was determined. 

In the canonical ensemble, several thermodynamic properties were investigated. Probability 
density functions for ^^•^''Fe were also extracted which reveal the difference between small 
systems such as the atomic nucleus and a system in the thermodynamic limit. 
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Microscopic model calculations based on BCS quasiparticles were performed. The overall 
agreement between the experimental and calculated level densities are good. Step structures 
observed in the experimental level densities are also observed in the plot of the average 
number of broken pairs as a function of excitation energy. The parity distributions obtained 
from model calculations for ^^'^^Fe indicate a decrease of parity asymmetry for both isotopes 
as the excitation energy increases. The microscopic model, which is expressed within the 
microcanonical ensemble (fixed E with no heat bath) , describes the observed level densities 
without the need of attenuation of the pairing-gap energies for excitation energies below 8 
MeV. 
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